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ABSTRACT 

A model for the chemical evolution of the intergalactic medium (IGM) is 
presented using theoretical yields of very massive (Mvms > 100 Mq) stars (VMSs) 
and Type II supernovae (SNe II). It is shown that if [Si/C] is indeed as high as 
~ 0.7 in the IGM, then VMSs (Mvms ~ 140-260 Mq) associated with pair- 
instability supernovae (PI-SNe) in low-mass (~ 10^ Mq) halos at high redshift 
must produce at least 50% of the Si. The remainder is from later galactic outflows 
of SN II debris, which also provide most of the C and O. Both sources are required 
to account for the metal inventory in the IGM. The early VMS production must 
continue until redshift z ~ 15 so that the efficiency of VMS formation per low- 
mass halo is significantly below unity. Contributions from the later galactic 
outfiows mainly occur at z ~ 4-6. Using a Salpeter initial mass function, we 
infer that the number of VMSs (Mvms ~ 260-2000 Mq) producing massive black 
holes (MBHs) with an average mass (Mmbh) ~ 270-550 is ~ 0.72 times the 
the number of VMSs associated with PI-SNe. The amount of metals (particularly 
Si) in the IGM that is attributable to PI-SNe is thus closely coupled with the 
total mass of MBHs produced in epochs prior to galaxy formation. Production 
of ~ 50% of the Si in the IGM by PI-SNe corresponds to an early inventory 
of MBHs that constitutes a fraction ~ (4-8) x 10~^ of the total baryonic mass 
in the universe. This is comparable to the global mass budget of the central 
supermassive black holes (SMBHs) in present-day galaxies. The corresponding 
occurrence rates in each halo of ~ 10^ Mq during the epoch of VMS formation 
at z > 15 are ~ 0.9 Gyr~^ for VMSs associated with PI-SNe and ~ 0.6 Gyr^^ for 
the concomitant more massive stars producing MBHs. These rates may be of use 
to studies of H2 dissociation and reionization and to models of SMBH formation. 

Subject headings: galaxies: formation — intergalactic medium — nuclear reac- 
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1. Introduction 

In this paper we explore implications of data on C, O, and Si in the intergalactic 
medium (IGM) for the "metal" sources. In earlier works (Wasserburg & Qian 2000; Qian 
& Wasserburg 2002; see also Qian, Sargent, & Wasserburg 2002), we proposed that very 
massive (Mvms > 100 M©) stars (VMSs) formed from big bang debris produced a prompt 
metal inventory in the IGM at the level of [Fe/H]iGM = log {Fe/U.)jQj^ — log (Fe/H)^ ~ —3, 
where (Fe/H) is the number ratio of Fe to H. This was based on observations of metal-poor 
stars in the Galactic halo that showed a sudden jump in the production of heavy rapid- 
neutron-capture (r-process) elements such as Ba and Eu relative to Fe at [Fe/H] ~ —3 (e.g., 
McWilliam et al. 1995; Burris et al. 2000). We inferred that the jump was due to the onset 
of major production of the heavy r-process elements by a class of supernovae associated with 
normal stars, which could predominantly form only after VMSs had provided a metallicity 
at the level of [Fe/H]iGM ~ —3 to the IGM. This is supported by theoretical models of star 
formation, which suggested that a critical metallicity of ~ 5 x 10~^ times the solar value is 
required for sufficient cooling of gas clouds to form low-mass protostellar aggregates (Bromm 
et al. 2001; see also Bromm & Loeb 2003). Using VMS yields of Heger & Woosley (2002, 
hereafter HW02), Oh et al. (2001) showed that there would be sufficient hard UV photons 
emitted from VMSs to reionize the universe if these stars were to provide [Si/HjjGM ~ —2.3 
(e.g., Cowie & Songaila 1998; Ellison et al. 2000; see also Aguirre et al. 2004). We here 
consider both early VMSs and later galactic outflows for metal enrichment of the IGM. It is 
argued that VMSs (Mvms ~ 140-260 Mq) associated with pair-instability supernovae (PI- 
SNe; HW02) are the dominant source of Si and also made substantial contributions to G 
and O. The dominant contributions to G and O are from later galactic outflows governed by 
Type II supernovae (SNe II). We then show that for a plausible initial mass function (IMF), 
the metal (particularly Si) inventory in the IGM due to PI-SNe is coupled with production 
of black holes by VMSs with Mvms > 260 M©. Thus, there is concomitant production of 
both metals and massive black holes at early epochs. 

The present study is an extension of our previous work (Qian & Wasserburg 2005, 
hereafter QW05), where we sought to explain the data on G, O, and Si in the IGM using a 
number of VMS and SN II models and considering a wide range of IMFs for VMSs. In that 
work, we showed from detailed arguments that VMSs associated with PI-SNe are required to 
account for the Si/C ratio in the IGM and that the preferred scenario for metal enrichment 
of the IGM is a combination of contributions from early VMSs and later galactic outflows. 
In the present paper, we assume the basic conclusions of that work and adopt the same input 
and formalism to develop the general self-consistent implications of the chemical evolution 
model used there. The new results to be presented include the constraint on the redshift for 
termination of VMS contributions and the relationship between metal production by PI-SNe 
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and the early inventory of massive black holes. 

We first review the issues. There is a substantial inventory of metals in the general IGM 
over a wide range of redshift, z — 1.5-5.5 (e.g., Songaila 2001; Pettini et al. 2003). Detailed 
studies of quasar absorption spectra have established abundances of the ionic species C III, 
C IV, O VI, Si III, and Si IV in the Lya forest (see Schaye et al. 2003; Aguirre et al. 2004; 
and Simcoe et al. 2004 for recent analyses). The corresponding elemental abundances have 
been determined using models of the UV background (UVB; Haardt & Madau 1996, 2001) 
to account for the fraction of an element in a specific ionization stage. Two typical UVB 
models are model Q, which has a harder UVB provided by quasars only, and model QG, 
which has a softer UVB provided by both quasars and galaxies. For both models, Schaye et 
al. (2003) and Aguirre et al. (2004) found that C and Si abundances are quite variable over 
different regions. For a given density of the Lya forest, the abundance of, say C, appears to 
follow a lognormal distribution. The net inventory of C over the density range (0.32 to 63 
times the cosmic mean) covered by the data is [C/H]igm = —2.3 for UVB model Q and —2.8 
for model QG. For either model, Schaye et al. (2003) and Aguirre et al. (2004) found that 
neither [C/H]igm nor [Si/H]iGM show evidence of evolution in the range z = 1.8-4.1. Using 
the data at 2; ~ 2.5 and similar UVB models, Simcoe et al. (2004) also found lognormal 
distributions for C and O abundances. The results of the above groups for [C/H]igm are in 
good agreement for either UVB model. The observed variability of abundances in different 
regions must be a result of varying degrees of mixing. However, there is no apparent variation 
of the Si/C or C/0 ratio (Aguirre et al. 2004; Simcoe et al. 2004). In the models presented 
below, we will focus on the net IGM inventory and the associated elemental ratios. The 
distribution of metals is not treated. 

From the observations there is no apparent evolution of the IGM inventory over z — 1.8- 
4.1. Thus, this inventory must have been produced prior to 2; ~ 4. The plausible sources 
are VMSs and SNe II. The Si/G and G/0 ratios^ are potential diagnostics of the sources. 

These ratios are sensitive to the UVB model. As both of these ratios are available for models 
Q and QG, we discuss these two UVB models first. While model Q gives [Si/C]iGM = 1-45 
and [C/0]iGM = 0, model QG gives much lower values of [Si/C]iGM = 0.74 and [G/0]igm = 
—0.50 (Aguirre et al. 2004; Simcoe et al. 2004). For both models [Si/G]iGM is quite 



^For both the IGM data and the results from stellar models, we use the same reference solar abundances 
as in QW05: loge0(Si) = log(Si/H)0 + 12 = 7.55 (Anders & Grevesse 1989) adopted by Aguirre et al. 
(2004), and loge0(C) = 8.52 and logeo(O) = 8.83 (Grcvcssc & Sauval 1998) adopted by Simcoe et al. 
(2004). As Schaye et al. (2003) and Aguirre et al. (2004) adopted a slightly different solar C abundance of 
loge0(C) = 8.55 (Anders & Grevesse 1989), we have made appropriate small adjustments to their [C/H]igm 
and [Si/C] IGM values for consistency. 
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high and points to VMS sources (Schaerer 2002; Aguirre et al. 2004). For UVB model 

QG, [Si/C]iGM is well above the yield ratio of SN II sources but below that of pure VMS 
sources (see Table 2 of QW05). The value of [Si/C]iGM for model Q is so high that it 
can be barely matched by the yield ratio of a pure VMS source. We have shown earlier 
(QW05) that the values of [Si/C]iGM and [C/0]igm for UVB model Q cannot be accounted 
for simultaneously by any existing stellar models. However, both elemental ratios for model 
QG can be matched using a combination of VMS and SN II sources. Results on [Si/C]iGM 
only were also given by Aguirre et al. (2004) for models QGS and QGS3.2, both of which 
have an even softer UVB than model QG. Model QGS gives [Si/C]iGM = 0.23 while model 
QGS3.2 gives [Si/C]iGM = 0.43. These [Si/C]iGM values are consistent with the yield ratio 
of SN II sources but significantly below that of pure VMS sources (see Table 2 of QW05). 
Consequently, were the composition of the IGM characterized by these [Si/C]iGM values, 
there would be no requirement of VMS contributions. However, Schaye et al. (2003) showed 
that for UVB model QGS, the C abundance strongly decreases with time for all densities, 
which is clearly unphysical. Furthermore, model QGS3.2 predicts a jump in the optical depth 
of Si IV relative to G IV due to the sudden transition in its UVB a,t z — 3.2. This is disfavored 
by the observations (Aguirre et al. 2004). In contrast, UVB model QG appears compatible 
with all the observations so far and also gives physical results regarding the evolution of 
metallicity with time and density (Schaye et al. 2003; Aguirre et al. 2004). Therefore, it 
seems unlikely that this model could be seriously in error (Aguirre et al. 2004). The results 
on the IGM inventory for UVB model QG (see Table 1) will be used in the subsequent 
presentation. As shown by QW05, VMS contributions are required to match [Si/C]iGM for 
this model and both [Si/C]iGM and [C/0]igm for this model can be accounted for by a blend 
of VMS and SN II contributions. Insofar as [Si/G]iGM is approximately at or above the 
value given by this model, these results appear robust. The high value of [Si/G]iGM was not 
discussed in the recent critique of VMS models by Tumhnson, Venkatesan, & ShuU (2004) or 
in another study by Daigne et al. (2004). In particular, Daigne et al. (2004) concluded that 
"there is no evidence, nor any need, for a hypothesized primordial population of very massive 
stars in order to account for the chemical abundances of extremely metal-poor halo stars 
or of the intergalactic medium." We consider that further improvements in measurements 
and analyses of the IGM abundances, as well as further theoretical studies of the conditions 
under which stars of extremely low metallicity can form, will clarify the choice between the 
metal enrichment scenarios with and without VMSs. 

Based on the high value of [Si/C]iGM for UVB model QG, we use the VMS model of 
HW02 and the low-metalhcity SN II model of Woosley & Weaver (1995, hereafter WW95) 
to quantify VMS contributions to the IGM. The VMSs relevant for metal enrichment are 
in the mass range Mvms ~ 140-260 and produce PI-SNe (HW02). A Salpeter IMF is 
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assumed for both VMSs and SNe II and the corresponding yield ratios are given in Table 1 
(models 2 A and 6 in Table 2 of QW05). Using these yield ratios, we can find mixtures of 
VMS and SN II contributions that can account for [Si/C]iGM and [C/0]igm simultaneously 
As both [C/0]vMS = -0.57 and [C/0]sn ii = -0.42 are within 0.1 dex of [C/0]igm = -0.50, 
the appropriate mixture is determined by matching [Si/CjicM. The required fraction /^^^ 
of Si contributed by VMSs can be calculated from 

V'^VIGM V'^VVMS V'^VSNII 

With [Si/C]vMS = 1-18 and [Si/C]sN ii = 0.42, ^ o.63 is required to match [Si/C]iGM = 
0.74 exactly and /^^^ = 0.5 gives [Si/C] = 0.65. Therefore, VMSs must be the dominant 

source of Si. This result does not depend on the specific stellar models adopted here. In 
fact, as other SN II models have lower values of [Si/C]sN ii (see Table 2 of QW05), using 
them only strengthens the above result. We conclude that a minimum value of f^j^^ ~ 0.5 
is required to account for the observed elemental ratios in the IGM. For the stellar models 
adopted here, /gX^^ _ o.5 corresponds to /^^^ ~ 0.15 and /^^^ ~ 0.19, so SNe II are the 
dominant source of C and O. This conclusion can be generalized to all stellar models as well. 



2. Metal Enrichment of the IGM by VMSs 

We now discuss metal enrichment of the IGM by VMSs. The formation of these stars 
is supposed to precede that of lower-mass stars and occur in halos where big bang debris 
can cool by II2 molecules (see Abel et al. 2002; Bromm & Larson 2004 for reviews). This 
requires a minimum halo mass of ~ 10^ Mq. The efficiency of VMS formation in such low- 
mass halos determines how early VMSs could provide substantial enrichment to the IGM. 
For a simple estimate, we consider that all halos with mass M > -Mmin are formed from 
"building-block" halos with mass Mmin. We assume that a VMS would have formed in every 
"building-block" halo and ejected all its nucleosynthetic products into the IGM. Then at 
redshift z, the abundance of element E in the IGM would be 

/7IGM _ (E/H)iGM ^ {Y^nnM > M^^\z) 
- (E/H)^ ~ XihM^^ ' 

where (V^^*^'^^) is the average mass yield of E per VMS event, F{M > M^i^\z) is the fraction 
of all matter residing in halos with M > Mmin at redshift z, is the solar mass fraction 
of E, and /& is the baryonic fraction of the halo mass. Consider Si as an example. Using 
(ygVMS^ ^ 16.6 Mo (HW02), X| = 7.1 x 10"^ (Anders & Grevesse 1989), fb = 0.15, and 
Mmin = 10^ Mo, we obtain Zlf^ ?s 1.6F(M > IO^MqI^). The function F{M > IO^MqI^) 
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can be estimated using the Press-Schechter formalism (Press & Schechter 1974) and ranges 
from ~ 0.3% to 3% for ^ ~ 24 to 17. Thus, if VMSs are to provide 50% of the measured 
inventory of [Si/H]iGM = log^gf^ = "2-0 (Aguirre et al. 2004), this requires that a VMS 
form in every "building- block" halo by ^ ~ 24 or more plausibly, in ~ 10% of such halos by 
17. 

For a formal and self-consistent model of the chemical evolution of the IGM using the 
rate of occurrence of a VMS event in a halo, we consider a large reference region of the 
universe and treat it as a closed homogeneous system (QW05). The evolution of (E/H)ig]v[ 
in the system as a function of time t is governed by 

c^(E/H)iGM ^ 

dt (H)iGM ' 

where Qe is the net rate for ejection of E atoms into the IGM and (H)igm is the number of 
H atoms in the IGM. The rate Qe depends on the rate of occurrence of a VMS event and 
the efficiency of gas expulsion by the VMS explosion. We consider that a VMS may form in 
a halo when the virial temperature Tyir of the gas reaches a minimum value of Tvir,o = 300 K 
for H2 coohng. For z :$> 1, 

T„,.2n(^U4^)"7— ) K. (4) 

V1.22/ \io^MqJ V 10 y 

where /x is the mean atomic weight and n = 1.22 or 0.6 for a neutral or ionized gas, respec- 
tively, with a primordial composition of H and He. We assume that all the gas would be 
expelled from the halo by the VMS explosion if the gravitational binding energy of the gas 
-E'bi.gas is less than the VMS explosion energy E'exp (e.g., Bromm, Yoshida, & Hernquist 2003) 
and retained for -Bbi.gas > -Sexp- For 2; » 1, 

E,,,^ ^ 4.31 X 10^^ (-^) ('4^) erg. (5) 



For a given z, a halo associated with an uq a density fluctuation (an no cr halo) would have 
the mass Mq corresponding to Tyir = Tvir,o for the onset of VMS formation, while an a 
halo would have the mass Mm corresponding to i?bi,gas = -Ecxp for the onset of gas retention. 
Then wc have Qe = Yln^no o-^E.n a-, where Nn a is the total number of halos with a 
specific n value in the system and Pe^u a is the number production rate of E in a single n a 
halo. 

The rate Pe^u a is related to the rate of occurrence of a VMS event i?„ „ in aiin a halo. 
We assume that these rates are proportional to the number (H)„ o- of H atoms in the gas of 
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the halo and write 

PE,na = (|/r')i?n . = Ar'(E/H)Q(H)„ (6) 
p ^ AVMSyO fbMn a /-v 

where (i/e'^^) is the average number yield of E per VMS event corresponding to the mass 
yield (1^^'^^) and A^^^ is a rate constant. The above prescription of P^^n a Rn a oc (H)„ g- 
is meant to imply that it is more probable for a VMS to form in a larger halo, and we will 
discuss the limit on the occurrence of a VMS event in a typical halo below. Substituting Qe 
in equation (3) and noting that for 2; ^ 1 the majority of the H atoms of the system resides 
in the IGM and the majority of the H atoms in a halo resides in the gas, we obtain 



= iVn .(H)„ ./(H)iGM « Ar'i^(Mo < M < Mbil^), (8) 



dt 

n=no 

Zf^{t) ^ A^^^is r^'^ p^M^ ^^j^ Mu\z')dt', (9) 







where 



F(Mo < M < Mu\z) = ^ J^^'^xp (^-^^ dx (10) 

is the fraction of all matter residing in n o" halos with no < n < nu as given by the Press- 
Schechter formalism and t{z) = 0.538[10/(1 + z)]^^^ Gyr for the adopted cosmology (see 
discussion in Barkana & Loeb 2001). 

The evolution of is thus only dependent on A^^'^ and F{Mo < M < Mu\z). As 

scales with A^'^^, we will use A^'^^ = 0.1 Gyr~^ as a reference value (for elements 
produced by SNe 11 only, the average Galactic value of Al'"^ is ~ 0.1 Gyr~^ over a period 
of ~ 10 Gyr prior to solar system formation; QW05). To determine F{Mq < M < M]^[\z), 
we use Tvir,o = 300 K (/x = 1.22) to find Mq and a typical VMS explosion energy of E^exp = 
4 X 10^^ erg (HW02) to find Mu for which Eu,gas = -E'exp- The relevant range of n cr halos 
with Mo < M < Mbi (no < n < nu) corresponds to 5.5 x 10^ Mq < M < 6.1 x 10^ 
(2.3 <n< 3.8) at z = 20 and 8.3 x 10^ Mq < M < 7.2 x 10^ Mq (1.8 < n< 2.9) at z = 15. 
Although the range Mq < M < is rather wide, the integral for F(Mo < M < Mki\z) is 
dominated by the contributions from n a halos with M ~ Mo (n ~ uq; see eq. [10]). Thus, 
a typical halo relevant for metal-enrichment of the IGM has M ~ 10^ M© for z ~ 15-20 and 
an exact treatment of halos of much higher masses, which requires a detailed description of 
their merger history, is not crucial to our model. [In particular, -F(Mo < M < Mbil-z) is not 
sensitive to ii^exp for E^^p > 10^^ erg as halos with -Ebigas ^ 10^^ erg are extremely rare for 



-8- 



^ > 1]. The evolution of [E/H]igm = log^^^^ for A^^^ ^ g.l Gyp-^ is shown as the dot- 
dashed curve in Figure 1. This gives [E/H]igm = —3.46 at z = 15. Thus, to provide 50% of 
the measured inventory of [Si/H]iGM = —2.0 by z = 15 (point A) requires Ag'j'^^ = 1.4 Gyr~-^, 
which corresponds to i?„ ^ti^ ~ 0.2 VMS event over ti^ = t{z = 15) = 0.27 Gyr for a halo 
with M„ o- ~ 10^ M0 and Rn a ^ 0.9 Gyr~^ (eq. [7]). In Figure 2 we show the values of 

^VMS 

necessary to achieve [Si/H] = —2.3 (solid curve) or —2.0 (dashed curve) by a given z 
as calculated from equation (9). The dot-dashed curve indicates formation of at{z) ~ 1 
VMS by time t{z) in a halo of mass M„ ~ 10^ Mq and is taken as a bound. It is evident 
that a major part (> 50%) of the Si in the IGM, which must come from VMSs, can only 
be provided for z substantially below 19 and with plausible production rates, for z < 16. 
We consider Ag^j^^ = 1.4 Gyr~^ a reasonable rate and show the corresponding evolution of 
[Si/H]iGM foi' -2 > 15 as the solid curve in Figure 1. The evolution of [C/H]igm and [0/H]igm 
is fixed by the yield ratios of the VMS model (see Table 1) and shown as the dot-dashed and 
short-dashed curves, respectively. 



3. Termination of VMS Contributions 

Figure 1 shows that for a reasonable formation efficiency, VMSs would provide [C/H] = 
—3.5 and [0/H] = —2.9 to the IGM hj z = 15. This is consistent with the proposal by 
Bromm & Loeb (2003) that VMSs can no longer be formed when [C/H] = —3.5 ± 0.1 and 
[0/H] = —3.05 ± 0.2 are reached in the IGM. They showed that for these metallicities, 
gas clouds initially cooled by H2 molecules could continue to cool by C and O atoms and 
fragment into smaller clumps. However, the termination of VMS formation is a complex 
matter. The soft UV radiation from VMSs would dissociate H2 molecules in low-mass halos, 
thus suppressing later VMS formation. The history of H2 dissociation is not well known 
but a theoretical study by Ciardi et al. (2000) suggested that universal dissociation has not 
yet occurred for z ~ 20. The UVB produced by VMSs may also play an important role 
in reionization, which could have taken place at 2; ~ 17 (Kogut et al. 2003). As shown 
by Oh et al. (2001), there should be sufficient hard UV photons emitted from VMSs to 
reionize the universe if these stars were to provide [Si/H]iGM ^ —2.3 (see Tumlinson et al. 
2004; Daigne et al. 2004 for a different view). The above estimates for the onset of C and 
O cooling, universal H2 dissociation, and reionization appear to roughly coincide. Clearly, 
there must be a transition region in which ongoing VMS formation increases the soft UVB 
for H2 dissociation and begins to suppress further VMS formation, thus resulting in a global 
decrease in Ag''^^. There will also be an accompanying increase in the hard UVB that 
begins to initiate reionization and raise the IGM temperature. The balance between these 
processes is not well understood and cannot be addressed here. It is possible that both effects 
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are operating to end the regime of H2 cooling at roughly similar z to reionization. Combining 
the consideration of the increasing UVB with the metallicity condition for termination of 
VMS formation proposed by Bromm & Loeb (2003), we will assume that VMS contributions 
to the IGM will cease at ^ ~ 15 and transition to formation of lower-mass stars will occur 
universally. 



4. Enrichment of the IGM by Galactic Outflows 

For simplicity we assume that VMS production stops sharply at 2; = 15 and take VMS 
contributions to the IGM as shown in Figure 1. Then the bulk of the C and O and the 
remainder of the Si must be provided by SNe 11 at 15 > 2; > 4. The formalism resulting 
in equation (8) can be modified to treat the chemical evolution of the IGM in this regime, 
where formation of regular stars including SN II progenitors requires Tvir,o = 10'^ K = 0.6) 
for coohng by atomic species. This corresponds to halos with M ~ 10^ Mq and -Ebi.gas far 
exceeding the typical SN II explosion energy of ~ 10^^ erg. Thus only a fraction e of the 
debris from an SN II will be ejected into the IGM. While the dependence of e on the halo 
mass may be complicated, we will treat it as a constant for halos below some cut-off mass 
Ml. Our goal is to estimate how efficient SN Il-driven galactic outflows must be in order 
to provide a substantial part of the IGM inventory at 15 > z > 4. Equation (8) may be 
rewritten in this regime as 

^ ^ eAf "F(Mo<M<Mi|^), (11) 



dt 



4«^(i) ^ Zf^{h,) + eKf r\{M^<M<M,\z')dt', (12) 

where ^^^^(^15) is the IGM inventory of E at 2; = 15 resulting from just VMS production. 

Note that Z^^{t) - ^^^^(^15) scales with eA|^ For convenience we use eA|^ " = 
0.01 Gyr~^ as a reference value. This corresponds to an outflow efficiency of e = 0.1 for 
a Galactic SN II production rate of " = 0.1 Gyr'^ Using = 10^ K and Mi = 
10^° M0, we show the contributions from just galactic outflows as the long-dashed curve in 
Figure 1. The value of Afi = 10^'^ Mq is consistent with observations of outflows from Lyman 
break and dwarf galaxies (e.g., Pettini et al. 2001; Martin et al. 2002). The results will not 
change much if Mi is varied from this value by a factor of several either way (QW05). As 
can be seen from Figure 1, galactic outflows can only provide [E/H]igm = —2.92 by 2; = 4 
for eA|N " = 0.01 Gyr-^. To provide [Si/H] = -2.3 (50% of the IGM inventory) by z = 4 
would require eAgf " = 0.042 Gyr~^. This corresponds to a rather high outflow efficiency 
of e = 0.42 for intermediate-mass (~ 10^-10^° Mq) halos with a Galactic production rate. 
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If galactic outflows are to provide [Si/H] = -2.3 by ^ = 3, then eA|f^ " = 0.024 Gyr"^ is 
required and e = 0.24 is sufficient for a Galactic production rate. 

Figure 1 shows an example of a self-consistent model for the chemical evolution of the 
IGM with contributions from VMSs for z > 15 and galactic outflows for z < 15. The overall 
evolution of [Si/H]iGM, [0/H]igm, and [C/H]igm is shown assuming that VMSs provided 

50% of the Si inventory hj z = 15 (point A) and galactic outflows the other half hj z = 4 
(point B) and using the yield ratios of the sources (see Table 1). It can be seen that galactic 
outflows cannot provide signiflcant contributions until 2; ~ 6. Thus, the model implies that 
there should be very little change in the net IGM inventory between z = 15 and 2; ~ 6. The 
respective contributions to C, O, Si, and Fe from VMSs and galactic outflows associated 
with SNe II are given in Table 1. 



5. Metals in the IGM and Inventory of Meissive Black Holes 

In considering metal enrichment of the IGM at high z we have focused on VMSs with 
^VMS ~ 140-260 Mq in low-mass halos using the yields of IIW02. These workers showed 
that this narrow range of stellar masses produce PI-SNe, which would efficiently eject the 
associated nucleosynthetic products into the IGM. They also showed that all VMSs with 
-^VMS > 260 Mq would form black holes and therefore not contribute any metals. Thus, the 
first stellar sources contributing metals are restricted to PI-SNe. However, if one considers 
any plausible IMF, the occurrence of PI-SNe implies the production of more massive stars 
that will thus give massive black holes (MBHs). It follows that the IGM inventory of metals 
that is attributable to PI-SNe must be coupled with the production of MBHs. This then 
relates the metal (particularly Si) content of the IGM to the inventory of MBHs at early 
epochs. 

Before discussing the coupling between the amount of metals in the IGM and the total 
mass of MBHs produced in epochs prior to galaxy formation, we note the work by Madau 
& Rees (2001), who estimated the inventory of MBHs based on considerations other than 
metal enrichment. Motivated by simulations of Bromm, Coppi, & Larson (1999) and Abel, 
Bryan, & Norman (2000), which suggested that the first stars may have been very massive, 
Madau & Rees (2001) recognized the significance of possible early production of MBHs in 
an episode of pregalactic star formation using earlier stellar models (Bond, Arnett, & Carr 
1984; Fryer, Woosley & Heger 2001), which showed that VMSs could produce MBHs. In 
particular, they found that if one MBH of mass AfMBH ^ 150 formed in each 3 a halo 
(with 3 X lO^h^^ Mq in dark and baryonic matter) collapsing at ^ ~ 20, then a fraction 
/mbh ^ 8 X 10~^h^ of all baryonic matter would be in MBHs, where h is the Hubble 
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parameter in units of 100 km Mpc~^ and we take h — 0.7 throughout this paper. This 
estimate is comparable to the total mass fraction of the central supermassive black holes 
(SMBHs) found in most nearby galaxies (see eq. [18] below). 

In the earher sections of the present paper, we have inferred that VMSs with Mvms ~ 
140-260 M© provided [Si/H]iGM = -2.3 (~ 50% of the Si in the IGM) and we have derived 
the rate of VMS formation required to achieve this. We now explore the implications of 
this metal production for the inventory of MBHs that would be produced in low-mass halos 
at the same epoch and show that it is in accord with the possibility considered by Madau 
& Rees (2001). For a quantitative estimate, we consider a large reference region of the 
IGM with a total mass Migm and use a Salpeter IMF for VMSs. To enrich this IGM with 
[Si/H]iGM = log Zlf^, the required number A/^vMS,Si of VMSs with Mvms ~ 140-260 JV/q is 

-vIGM M 

^VvMS.Si ^yVMS^ ■ V^-^j 

As the total baryonic mass of a hosting halo is ~ 10^ Mq, we take the upper mass hmit for 
VMSs to be 2000 Mq. We have checked that our results are not sensitive to this assumed 
upper mass limit. The number A^mbh of MBHs resulting from VMSs with Mvms — 260- 
2000 Mq is 

NuBK ~ ^tIo TTT— ^VMS,Si ~ 0.727VvMS,Si, (14) 

ji40 dm 

where m is the VMS mass in units of Mq. The average mass of these MBHs is 



(Mmbh) ~ «^S)o ^^T^^o ~ 547aM0, (15) 



where a ~ 0.5-1 is the average mass fraction of a VMS progenitor in the relevant mass range 
that ends up in the MBH (Fryer et al. 2001). The fractional contribution of MBHs to the 
baryonic mass of the universe can be estimated as 

„ _ A^MBh(Mmbh) _ ^IGM^0 [ ^MBH ^ ( (Mmbh) ^ 

The above estimate^ of /mbh may be compared with the fraction /smbh of all baryonic 



^As stated above, this result assumes a Salpeter IMF (cx m ^■^^) for VMSs. For a steeper IMF of the 
form cx m~^, the numerical coefficient in eq. (17) for /mbh is reduced by a factor of w 2. 
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matter contributed by SMBHs in the present universe, 

/SMBH = X 10-^ (18) 

^ ''bp cxi 

where psmbh = 4.6li;4 ^ 10^ M© Mpc~^ is the present mass density of SMBHs obtained 

from galactic observations (sec the recent analyses by Marconi et al. 2004), pcri = 1-36 x 
10"*^^ Mq Mpc^'^ is the present critical density, and Vti, = 0.045 is the fractional contribution 
of baryons to p^ri- It can be seen that the mass fraction of MBHs calculated from the Si 
abundance in the IGM of the early universe (eq. [17]) and that of SMBHs determined for 
the present epoch by observations (eq. [18]) are in remarkably good accord. The possibility 
considered by Madau & Rees (2001) is now shown to be supported by the required VMS 
contribution to produce the Si in the IGM. We note that Schneider et al. (2002) derived 
a wide range of formation efficiency for VMSs with Mvms ~ 140-260 Mq using /mbh < 
/sMBH and other constraints. Their results can be easily accommodated by the rate of 
VMS formation derived here. We also note that accreting MBHs may contribute to the 
soft X-ray background (SXRB). Based on the unaccounted SXRB flux, Salvaterra, Haardt, 
& Ferrara (2005) estimated that /mbh cannot exceed ~ 10^'^. This is consistent with our 
result. During the epoch of VMS formation at 2; > 15, the rate of occurrence of an MBH- 
producing VMS event is (-/VMBH/-^VMS,Si)-Rn a ~ 0.6 Gyr"-*^ for a halo of ~ 10^ (eq. [7] 
with Ag^^^ = 1.4 Gyr~^ and eq. [14]). This rate can be used to calculate the contribution 
from MBHs to the SXRB. 

The remarkable accord between the MBH inventory calculated from the Si content of the 
early IGM (eq. [17]) and the present global mass budget of SMBHs obtained from galactic 
observations (eq. 18]) cannot be easily dismissed as an accident. It is expected that the 
early-formed MBHs would cluster near galactic centers (Madau & Rees 2001). As a fraction 
~ (4-8) X 10~^ of all baryonic matter is in MBHs (eq. [17] with a ~ 0.5-1), a typical 
galaxy with a baryonic mass of ~ 10^^ Mq would have a total mass of ~ (4-8) x 10^ Mq 
in MBHs. This is close to the mass of (3.7 ± 0.4) x 10^ M© for the SMBH at the Galactic 
center (e.g., Ghez 2004). It is possible that the central SMBHs in typical galaxies of the 
present universe are formed by simply merging the MBH seeds without much gas accretion. 
This would explain why the total mass of SMBHs observed at the present epoch corresponds 
approximately (within a factor of ~ 2) to that of the MBH seeds. 

However, there appears to be an obvious conflict. It is widely accepted that quasars are 
powered by gas accretion onto SMBHs. The total gas mass accreted by all optical quasars 
is inferred to be also comparable to that of the SMBHs observed at the present epoch 
(e.g., Yu & Tremaine 2002; Marconi et al. 2004). This is supported by the simulations 
of Hopkins et al. (2005), who proposed a unified model for producing SMBHs, quasars, 
and galaxy spheroids. In their model, SMBHs were formed predominantly by gas accretion. 
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Both the global mass budget and the relative distribution of SMBHs over the mass range 
-^SMBH ~ 10^-10^° Mq given by their model are in agreement with those derived by Marconi 
et al. (2004) from galactic observations. It appears that the model of Hopkins et al. (2005) 
is rather satisfactory. However, within the uncertainties of their model, it is still plausible to 
consider that the more common SMBHs of < 10^ Mq might represent aggregates of the MBH 
seeds without much gas accretion but large amounts of gas accretion must occur to explain 
the rarer SMBHs of > 10^ Mq. Thus, in accounting for the approximate agreement between 
the global SMBH mass budget determined from galactic and quasar observations and the 
MBH inventory, there appears to be a "rule" that only the rarer higher-mass SMBHs can 
efficiently accrete gas but the more common lower- mass ones cannot. This matter clearly 
requires attention. In the model of Hopkins et al. (2005), mergers of halos enhance the infall 
of gas toward the SMBHs. However, as mergers must have also occurred to form galaxies 
hosting SMBHs of < 10^ M©, which constitute ^ 50% of the global SMBH mass budget 
(Marconi et al. 2004), there seem to be some other processes that inhibit gas accretion onto 
such lower-mass SMBHs. 

The complex problem of forming SMBHs from MBH seeds was studied by Volonteri, 
Haardt, & Madau (2003). However, the total mass of MBHs in their scenarios is only a 
small fraction (typically ~ 10~^) of that of the SMBHs at the present epoch as cited above. 
We urge that the problem of forming SMBHs from MBH seeds be revisited with the initial 
conditions presented here. For a formation rate of ~ 0.6 Gyr"^ in each halo of ~ 10^ M©, 
an MBH with an average mass (Mmbh) ~ 270-550 Mq (eq. [15] with a ~ 0.5-1) should 
have been produced in ~ 16% of such halos by 2; = 15. If the amount of gas accretion 
onto black holes in galactic centers is typically comparable to or smaller than the total mass 
of the original inherited MBHs, then there would not be a serious difficulty in reconciling 
the global SMBH mass budget determined from galactic observations with that from quasar 
observations and with the inventory of MBHs calculated from the VMS contribution to the 
Si in the IGM. 

The evolution of an initial distribution of MBH seeds was studied extensively by Islam, 
Taylor, & Silk (2003, 2004) using a semi-analytical model to follow the hierarchical merging 
of halos. These authors focused on the dynamics of the MBH inventory and presented a 
clear analysis of what might be expected from assemblage of MBHs without gas accretion. 
Amongst their conclusions are: (1) Hierarchical merging of MBH seeds that formed in ~ 3 cr 
halos collapsing at 2; ~ 25 can contribute > 10% of the present global mass budget of 
SMBHs. A central SMBH of ~ 3 x 10^ Mq in a galaxy hke ours could be the result of 
seed accumulation without gas accretion. However, gas accretion is necessary for higher- 
mass SMBHs. (2) For a present-day galaxy of ~ lO^'^-lO'^^ Mq (mostly in dark matter), 
the total mass of the inherited MBHs in the galactic halo is comparable to or greater than 
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the mass contributed by MBHs to the central SMBH. It would be of interest to see what 
a new calculation of the type presented by Hopkins et al. (2005) would give if the mass 
and production history of MBHs as provided here were used. If the majority of MBHs are 
stored in galactic halos but not bulges, then we might expect gas accretion to account for 
the growth of central SMBHs while the MBHs in the low-density halos might not accrete. In 
this case, it is accidental that the total mass of MBHs calculated from the VMS contribution 
to the Si in the IGM is close to the present global mass budget of SMBHs. 

6. Conclusions 

A two-stage model for the chemical evolution of the IGM is proposed considering VMSs 
in low-mass (~ 10^ M©) halos during early epochs {z > 15). This early stage ends when 
universal H2 dissociation and reionization occur and is followed by an extended quiescent 
period of little metal production until sufficient SNe II could occur in intermediate-mass 
(~ 10^-10^*^ Mq) halos to drive significant outflows. Galactic outflows contribute mainly 
at 2; ~ 4-6. The requirement of early VMS contributions follows from the high value of 
[Si/C] ~ 0.7 inferred for the IGM (Aguirre et al. 2004; see discussion therein regarding 
possible uncertainties) and model yields of VMSs and SNe II (e.g., HW02; WW95). The 
bulk (> 50%) of the Si must come from VMSs, while the C and O are predominantly from 
the later galactic outflows. The required VMS formation efficiency corresponds to ~ 0.2 
VMS per low-mass halo. The requirement on galactic outflows implies efficient (~ 40%) 
loss of SN II debris from intermediate-mass halos with a Galactic SN II rate. Contributions 
from galactic outflows would be diminished as the halo masses and the amount of baryonic 
matter stored in low-mass stars increase. This may explain the lack of evolution of the IGM 
inventory for z — 1.8-4.1. 

The VMS contributions are consistent with the level of metaUicities at which VMS 
formation would be suppressed (Bromm & Loeb 2003). They also correspond to a sufficient 
number of UV photons from VMSs to reionize the universe (Oh et al. 2001). However, how 
the increasing photon production by VMSs leads to H2 dissociation and reionization, which 
then terminates VMS formation, is a complex issue and not addressable here. There is a 
hint that all these events occur at similar z. It would be a useful test if H2 dissociation and 
reionization could be modeled with the VMS formation rates proposed here. 

In producing [Si/H] = -2.3, VMSs also provide [Fe/H] = -2.9 to the IGM (see Table 1). 
This Fc abundance was considered by us (e.g., Wasscrburg & Qian 2000; Qian & Wasserburg 
2002) as the onset of regular star formation based on observations of low-metallicity Galactic 
halo stars. There is a basic issue regarding the metallicity at which regular stars may be 
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formed. As C and O atoms provide the important cooling, Bromm & Loeb (2003) proposed 
conditions for regular star formation in terms of threshold C and O abundances. These 
abundances coincide with the VMS contributions in our model. In this sense the threshold 
Fe abundance proposed by us is equivalent to the threshold C and O abundances. However, 
the low-mass star HE 0107-5240 cannot be explained by our scenario. This special star has 
[C/H] well above and [0/H] at or below the threshold but [Fe/H] = —5.3 (Christlieb et al. 
2004; Bessell, Christlieb, & Gustafsson 2004). Whether some low-mass stars such as this 
may form along with VMSs remains an open issue. 

A number of Galactic halo stars arc observed to have —4 < [Fc/H] < —3. If they formed 
at 2; < 4, this would pose a problem as our model gives a high value of [Fe/H]iGM = —2.3 
(see Table 1). However, the Fe abundances of these stars may be explained if they formed 
at 6 < 2; < 15. The IGM then would have a net inventory of [Fe/H] ~ —2.9 but the 
Fe abundances in different regions would follow a lognormal distribution with a scatter of 
~ 0.75 dex (Schaye et al. 2003; Simcoe et al. 2004). This leads us to suggest that some 
parts of the Galaxy formed from under-enriched regions of the IGM. 

In conjunction with the required metal production by VMSs in the narrow mass range 
-^VMS ~ 140-260 M0, there would also be MBHs produced by more massive stars (HW02). 
Both of these processes should occur in low-mass halos for z > 15. This leads to a relationship 
between the metal (particularly Si) content of the early IGM and the inventory of MBHs 
((Mmbh) ~ 270-550 Mq) at the same epoch. The total mass of MBHs calculated from 
the VMS contribution to the Si in the early IGM using a Salpctcr IMF is in remarkable 
agreement with that of SMBHs observed at the present epoch. These early-formed MBHs 
could cluster near galactic centers during the later epochs of galaxy formation when galactic 
outflows contributed to the IGM. Such a cluster of MBHs would, if coalesced, form an SMBH 
of ~ (4-8) X 10^ Mq in a typical galaxy with a baryonic mass of ~ 10^-*^ Mq. There is a hint 
of quasi-conservation of black hole masses, which implies that only the rarer SMBHs of 
> 10^ Mq can efficiently accrete gas but the more common ones of < 10^ Mq cannot. It 
remains to be investigated whether a mechanism exists to inhibit gas accretion onto lower- 
mass SMBHs. In conclusion, early "chemistry" appears to provide considerable insights 
into aspects of cosmological problems and offer many intriguing possibilities regarding larger 
cosmological issues. 
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Fig. 1. — Model for evolution of abundances of Si (solid curve), O (short-dashed curve), 
and C (dot-dashed curve). The sources are taken to be VMSs (yields from HW02) in low- 
mass halos ioi z > 15 and galactic outflows of SN II products (yields from WW95) from 
intermediate-mass halos for 2; < 15. The evolution is governed by the rates A^^^^ and eA|^ 
For Si these are chosen to give [Si/H] = —2.3 at 2; = 15 (point A) and the full IGM inventory 
of [Si/H] = —2.0 at 2; = 4 (point B). The rates for O and C are fixed by the yields of the 
sources. The long-dashed curve represents the evolution resulting from just galactic outflows. 
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Fig. 2. — Production rates of Si required to provide 50% (solid curve) and all (dashed curve) 
of the IGM inventory of [Si/H] = —2.0 by a given z as calculated from equation (9). The 
dot-dashed curve corresponds to the condition that every halo of mass M„ ^ ~ 10^ Mq would 
have Rn at{^) ^ 1 VMS by a given z. If VMS formation occurs at ~ 10% of the time in such 
hales, the bulk of the Si inventory can only be obtained at 2; < 15. 
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Table 1. IGM Data, Stellar Models, and Example Results^ 





[Si/C] 


[C/0] 


[0/Fe] 


[Si/H] 


[C/H] 


[0/H] 


[Fe/H] 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


IGM 


0.74 


-0.50 




-2.0 


-2.8 


-2.3 




VMSs 


1.18 


-0.57 


-0.02 


-2.3 


-3.5 


-2.9 


-2.9 


SNc II 


0.42 


-0.42 


0.17 


-2.3 


-2.7 


-2.3 


-2.5 


Mixture 


0.65 


-0.45 


0.13 


-2.0 


-2.7 


-2.2 


-2.3 



'^Data on the IGM inventory are for UVB model QG and are from 
Schaye et al. 2003, Aguirre et al. 2004, and Simcoe et al. 2004. Yield 
ratios for the VMS model of HW02 and the SN II model of WW95 
assuming a Salpeter IMF (models 2A and 6 in Table 2 of QW05) are 
given in cols. (2)-(4). The reference solar abundances are the same as 
adopted in QW05. A mixture is calculated where VMSs and SNe II 
contribute equally to the Si. The respective contributions from VMSs 
and SNe II are given in cols. (5)-(8). 



